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An  integrated  guidance  and  control  (IGC)  design  approach  is  proposed  based  on  small-gain  theorem  for 
missiles  steered  by  both  canard  and  tail  controls.  The  angle  of  attack  and  pitch  rate  commands,  which 
are  aimed  at  producing  desired  aerodynamic  lift  to  achieve  robust  tracking  of  a  maneuvering  target,  are 
generated  by  a  guidance  law  that  is  designed  using  input-to-state  stability  (ISS)  theory.  An  IGC  law  is 
developed  utilizing  generalized  small-gain  theorem  to  enforce  the  commands,  and  it  can  be  shown  that 
both  the  line-of-sight  (LOS)  rate  and  the  tracking  error  are  input-to-state  practically  stable  (ISpS)  with 
respect  to  target  maneuvers  and  missile  model  uncertainties.  The  algorithm  is  tested  using  computer 
simulations  against  a  maneuvering  target. 

©  2012  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

The  guidance  and  control  systems  of  interceptor  missiles  are 
usually  designed  separately  and  then  integrated,  but  it  can  be 
argued  that  this  scheme  cannot  fully  exploit  synergistic  relation¬ 
ships  between  the  two  subsystems  or  strictly  maintain  the  stability 
of  the  overall  system.  In  particular,  as  the  relative  range  between 
interceptor  and  target  becomes  small,  rapid  changes  in  the  rel¬ 
ative  geometry  might  cause  the  separation  design  method  to  be 
invalid  (Shima,  Idan,  &  Golan,  2006).  Integrated  guidance  and  con¬ 
trol  (IGC)  design  is  one  of  emerging  trends  in  missile  control  tech¬ 
nology,  because  it  can  improve  the  performance  of  interceptors  by 
viewing  guidance  and  control  loops  as  an  integrated  system  and 
taking  couplings  between  subsystems  into  account.  Such  a  design 
can  reduce  the  cost  of  the  required  sensors  and  increase  the  system 
reliability  (Williams,  Richman,  &  Friedland,  1983).  Since  the  canard 
fins  can  generate  an  aerodynamic  force  to  drive  the  missiles  to  the 
required  direction,  a  canard  configuration  was  employed  in  some 
short-range  air-to-air  missiles  (Shima  et  al.,  2006).  However,  ca¬ 
nard  control  has  its  limit  due  to  aerodynamic  saturation  at  high 
angles  of  attack,  especially  for  long  fuselage  missiles,  and  a  tail 
control  is  often  preferred.  Therefore,  some  researchers  proposed 
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IGC  approaches  by  using  both  canard  and  tail  controls  to  obtain 
better  performance  for  dual-control  missiles  (Idan,  Shima,  8t  Golan, 
2007). 

After  IGC  design  was  put  forward  in  Williams  et  al.  (1983), 
various  control  methods  have  been  introduced.  Suboptimal 
control  methods  were  introduced  to  solve  IGC  problem  in  three 
dimensions  (Menon  &  Ohlmeyer,  1999;  Palumbo  &  Jackson, 
1999;  Xin,  Balakrishnan,  8t  Ohlmeyer,  2006).  The  state  dependent 
Riccati  equation  (SDRE)  technique  was  employed  in  Menon  and 
Ohlmeyer  (1999)  based  on  a  six-degrees-of-freedom  (6-DOF) 
nonlinear  missile  model  for  IGC  design  with  the  assumption  that 
the  target  velocity  can  be  negligible  when  compared  with  the 
missile  velocity.  Palumbo  and  Jackson  (1999)  formulated  the  IGC 
problem  as  a  single  minmax  optimization  problem  when  the 
target  acceleration  was  known  as  a  priori  to  find  a  controller  that 
minimizes  the  final  miss  distance  and  control  energy,  in  which  the 
SDRE  technique  was  employed  to  handle  this  finite-time  horizon 
nonlinear  problem.  However,  as  mentioned  in  that  paper,  since 
solving  SDRE  online  is  time  consuming,  it  may  be  not  feasible 
for  a  6-DOF  missile  system.  Xin  and  co-workers  (Xin  et  al,  2006) 
solved  the  nonlinear  infinite-horizon  IGC  problem  by  utilizing 
0  —  D  technique.  Their  method  gave  an  approximate  closed-form 
suboptimal  feedback  controller  with  no  iterative  solutions  as  in  the 
case  of  the  SDRE  approach,  but  the  disturbances  and  uncertainties 
of  the  6-DOF  missile  model  were  not  considered. 

Sliding-mode  control  (SMC)  is  another  typical  method  that 
is  utilized  to  design  IGC  algorithm  in  pitch  channel  (Idan  et  al., 
2007;  Shima  et  al.,  2006;  Shtessel,  Shkolnikov,  &  Levant,  2009; 
Shtessel  &  Tournes,  2009),  especially  for  dual-control  missiles 
(Idan  et  al.,  2007;  Shtessel  &  Tournes,  2009).  In  Shtessel  et  al. 
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(2009),  a  sliding  surface  that  depends  on  the  line-of-sight  (LOS) 
rate  was  defined  in  the  guidance  loop  with  the  missile  pitch 
rate  viewed  as  a  virtual  control.  In  the  control  loop,  the  second- 
order  SMC  was  used  to  control  the  pitch  rate  to  track  the  virtual 
control  robustly  in  finite  time.  Shima  and  co-workers  (Shima  et  al„ 
2006)  used  SMC  to  obtain  IGC  design  by  using  zero-effort  miss 
(ZEM)  distance  as  a  single  sliding  surface  with  the  assumption 
that  the  target  acceleration  can  be  measured.  For  missiles  with 
both  canard  and  tail  controls,  they  chose  an  additional  sliding 
surface  based  on  flight-control  considerations  that  can  achieve 
improved  stability  and  shaped/damped  response  (Idan  et  al.,  2007). 
Note  that,  in  order  to  remove  nonlinear  terms,  the  equations  of 
IGC  model  in  Idan  et  al.  (2007)  and  Shima  et  al.  (2006)  were 
formulated  under  the  assumption  that  the  angle  between  LOS 
and  missile  velocity  is  almost  constant,  but  this  assumption  might 
be  not  proper  in  practice  since  large  maneuvers  of  a  target  may 
lead  to  significant  variation  of  that  angle.  In  Shtessel  and  Tournes 
(2009),  an  IGC  algorithm,  integrated  with  the  smooth  second-order 
sliding  mode  guidance  law  in  Shtessel,  Shkolnikov,  and  Levant 
(2007),  was  developed  using  higher-order  SMC  for  interceptors 
steered  by  a  combination  of  aerodynamic  lift,  sustainer  thrust, 
and  center-of-gravity  divert  thrusters.  With  this  approach,  a  fast 
tracking  of  the  attitude  command  can  be  achieved,  however,  it 
may  not  strictly  maintain  the  stability  of  the  overall  system  during 
the  process  of  tracking  commands  in  each  channel.  Backstepping 
scheme  was  also  used  to  design  IGC  algorithm  in  the  pitch  channel 
in  Sharma  and  Richards  (2004)  by  linearizing  the  missile  dynamics 
and  neglecting  the  lift  contribution  of  control  surfaces  based  on 
the  assumption  that  the  neural  network  can  provide  accurate 
estimates  of  uncertainties. 

Besides  the  works  mentioned  above,  the  feedback  linearization 
method  (Menon  &  Ohlmeyer,  2001),  subspace  stabilization  strat¬ 
egy  (Tournes  &  Wilkerson,  2001)  and  some  other  control  meth¬ 
ods  were  also  utilized  in  the  IGC  design  problem.  These  researches 
made  great  contributions  to  the  development  of  IGC  design,  but 
many  existing  results  were  obtained  based  on  some  strong  as¬ 
sumptions  or  without  considering  robustness  against  uncertainties 
and  disturbances. 

In  this  paper,  a  novel  IGC  design  approach,  against  maneuvering 
targets,  is  proposed  using  small-gain  theorem  (Jiang,  Teel,  & 
Praly,  1994)  and  input-to-state  stability  (ISS)  (Sontag,  1989)  for 
missiles  controlled  by  forward  and  aft  control  surfaces.  An  ISS- 
based  guidance  law  is  designed  for  target  interception  meaning 
that  the  law  keeps  the  LOS  rate  within  a  small  neighborhood  of 
zero  in  the  presence  of  unknown  target  maneuvers.  Then,  an  IGC 
law  is  developed  utilizing  small-gain  theorem  to  control  angle  of 
attack  and  pitch  rate  to  track  commands  generated  by  an  inversion 
depending  on  the  ISS-based  guidance  law.  Theoretical  analysis 
shows  that  the  IGC  approach  makes  both  the  LOS  rate  and  the 
tracking  error  of  attitude  angle  (rate)  be  input-to-state  practically 
stable  (ISpS)  with  respect  to  target  maneuvers  and  missile  model 
uncertainties  which  are  not  needed  to  be  known  as  a  priori.  It 
is  worth  claiming  that  our  approach  is  formulated  without  the 
assumption  that  the  angle  between  LOS  and  missile  velocity  is 
almost  constant,  and  the  lift  contribution  of  control  surfaces  is  also 
taken  into  account.  Besides  that,  the  stability  of  the  overall  system 
can  be  guaranteed  by  small-gain  theorem. 


Fig.  1.  Planar  interception  geometry. 


where  R  is  relative  range;  q  is  line-of-sight  angle;  <Pm  and  cpj  are 
missile  and  target  fight  path  angles,  respectively;  and  VM  and  VT  are 
missile  and  target  velocity,  respectively.  Suppose  that  the  missile 
speed  is  constant  during  the  end  game,  and  the  missile  acceleration 
can  be  calculated  as 

aM  =  VMipM.  (2) 

Differentiating  equation  (lb),  we  have 
_Vr  .  %  .  .Am 

0>q  =  -2— (0q+  -j±  -cos(q-<pM)—  (3) 

where  14  =  R,coq  =  q,  and  aTq  is  the  projection  of  bounded  target 
acceleration  orthogonal  to  the  LOS.  It  is  rational  to  assume  that 
VM  >  VT  and  \q  —  <pM\  <  |  when  homing  guidance  begins,  and 
because  a  guidance  law  should  nullify  the  LOS  rate,  \q  —  cpM\  <  § 
will  hold  generally  in  the  whole  process  of  homing  guidance,  that 
is,  there  exists  a  positive  real  number  of  k  such  that  cos (q  —  cpM) 


2.2.  Missile  dynamics 


Fig.  1  gives  the  missile  planar  dynamics,  where  C  is  the  missile’s 
center  of  gravity,  Cxb  is  aligned  with  the  missile’s  longitudinal  axis, 
and  Cxg  is  parallel  to  the  horizontal.  The  planar  missile  dynamics 
with  both  canard  and  tail  controls  are  given  by  Idan  et  al.  (2007) 

]io  =  M(a ,  co,  8C,  St)  (4a) 

b=co  (4b) 


I(or,  Sc,  St) 

<f>M  = 


(4c) 

(4d) 


where  J  is  moment  of  inertia;  co  is  pitch  rate;  &  is  pitch  angle;  m 
is  missile  mass;  a  is  angle  of  attack;  Sc  and  <5t  are  canard  and  tail 
controls,  respectively;  and  L  and  M  are,  respectively,  aerodynamic 
forces  and  moments,  which  are  nonlinear  functions  of  the  related 
variables,  in  particular  a,  co,  Sc,  and  8t.  It  is  assumed  that  the 
missile  actuators  follow  first-order  dynamics  with  time  constants 
tc  and  rt  as  expressed  by  Eq.  (4e).  When  a,  co,  Sc,  and  8t  are  in 
relatively  small  domains  around  zero,  Eqs.  (4)  can  be  transformed 
into  a  linearized  short-period  form  (Idan  et  al.,  2007;  Sharma  & 
Richards,  2004;  Shima  et  al.,  2006;  Wise  &  Broy,  1998),  that  is. 


la_ 

"  vMa 


(5a) 


2.  Model  derivation 

2. 1 .  Engagement  kinematics 

The  corresponding  equations  of  motion  between  missile  and 
target  are  given  by  Zhou,  Sun,  and  Teo  (2009) 

R  =  VT  cos (q  -  cpr)  -  VM  cos (q  -  cpM)  (la) 

Rq  =  -VT  sin (q  -(pT)  +  VM  sin (q  -  cpM)  ( lb) 


b  =co 

(5b) 

co  =  maOL  +  mwa 

)  +  mScSc  +  mSt8t  +  Am 

(5c) 

oc  =  ■&  -cpM 

(5d) 

8  S°c  ~  Sc 

8  “  St 

(5e) 

with  lx  =  3(^m) ,  mx  =  a(^J)  ,  and  Aa  and  Am  are  bounded  uncer¬ 
tainties. 
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3.  Integrated  guidance  and  control  law  design 


Applying 


3.1.  Concepts  and  preliminaries 

Consider  the  following  general  interconnected  system 
Hi  :  *1  =/i(xi,y2,Ui),  =  h1(x1,y2,ul)  (6) 

H2  :  x2  =  /2(x2,yi,u2),  y2  =  h2(x2,yuu2)  (7) 

where,  for  i  =  1,  2,  e  R"',  u,-  e  Rmi,  and y,-  e  RPi.  The  functions 
fi,f2,  hi  and  h2  are  smooth  and  a  smooth  function  h  exists  such 
that 

(yi,y2)  =  h(xi,x2,Ui,u2) 

is  the  unique  solution  of 

jyi  =  h,(xu  h2(x2,yu  u2),  ut) 

\y2  =  h2(x2,  hi(xi,y2,  Ui),  u2). 

We  have 


05i<  i  2 , 

q  R  ~  2|2  q  2  R2 
where  £  >  0,  into  Eq.  (12),  we  obtain 

+  ^  -«(1-  »)f) 

Substituting  (10)  into  Eq.  ( 14)  yields 


2  R2 


(13) 


(14) 


%2  aT0 

V  <-Kco2q  +  j-^.  (15) 

In  the  interception  endgame,  inequality  0  <  Rm  <  R  <  RM  holds. 
Hence, 


2  I2  ai 


(16) 


Theorem  1  (/iang  et  al„  1994).  Suppose  (6)  and  (7)  are  ISS  with 
(y2,  Ui)  (respectively  (yi ,  u2) )  as  input,  yi  (respectively  y2)  as 
output,  and  there  exist  class  X  £  functions  j3i,  f)2,  class  X  functions 
Vi’  Yi'  Yi'  Y2<  and  nonnegative  constants  d\,d2  such  that 

l|yi(OII  <  ft (.11*1  (0)||,  t)  +  yf(||y2t!l)  +  y“(||ui||)  +  di  (8a) 


Solving  the  differential  inequality  (16)  yields 

fc2 

V(oM(t))  <  e-2,aV(ojqm  +  ^  (1  -  e~2Kt) 

x  sup  a2q( r).  (17) 


lly2(t)||  <ft(||x2(0)||,t)  +  y2y(||ylt||)  +  y2u(||a2||)  +  d2.  (8b) 

If  two  class  Xoo  functions  p\  and  p2  and  a  nonnegative  real  number 
si  satisfying 


|  (Id  +  pi)  o  y2  °  0d  +  Pi)  o  y\ (s)  <  s 
{ (Id  +  pf)  o  yf  o  (Id  +  pi)  o  yy2  (s)  <  s  ’ 


Vs  >  Si 


(9) 


exist,  where  Id  denotes  the  identity  function  and  o  is  composition  of 
functions,  system  (6)-(7)  with  u  =  (ui ,  u2 )  as  input, y  =  (yi, y2)  as 
output  and  x  =  (xi,  x2)  as  state  will  be  input-to-output  practically 
stable  (IOpS)  (input-to-output  stable  (IOS)  if  st  =  di  =  d2  =  0J. 


3.2.  ISS-based  guidance  law  design 


Taking  the  square  roots  and  using  the  inequality  Va2  +  b2  <  a  +  b 
for  nonnegative  numbers  a  and  b,  we  can  see  that  Eq.  (11)  holds. 
Therefore,  the  closed-loop  system  of  system  (3)  and  guidance  law 
(10)  is  ISS  with  respect  to  the  target  maneuvers. 

Moreover,  if  the  target  normal  acceleration  vanishes,  that  is, 
aTq  =0,  Eq.  (11)  can  be  rewritten  as  \ooq(t)\  <  e^Kt\ojq(0)\.  In 
this  case,  the  origin  of  the  closed-loop  system  is  exponentially 
stable.  □ 

Remark  1.  In  guidance  law  (10),  —2VRa>q  stands  for  realistic  true 
proportional  navigation  (RTPN)  guidance  (Yang  &  Yang,  1997).  That 
is,  our  guidance  law  is  composed  of  a  RTPN  term  and  the  term  Rojq 
with  guidance  coefficients  K  and  £. 


Accepting  the  intuition  that  zeroing  the  LOS  rate  will  lead  to 
interception,  we  will  design  a  guidance  law  to  keep  the  LOS  rate 
within  a  small  neighborhood  of  zero  in  this  section. 

Theorem  2.  In  the  interception  of  a  maneuvering  target,  assume  that 
the  inequality  cos(q  —  <pM)  >  k  (k  is  a  positive  real  number) 
holds  in  a  reasonable  flight  domain.  Then  the  closed-loop  system  of 
system  (3)  and  guidance  law 

=  (,0) 

is  ISS  with  respect  to  the  target  maneuvers  for  K  >  0  and  f  >  0,  that 
is, 

\coq(t)\  <  e-KV,(0)|  +  — J^Vl-e-2®  sup  |aT,(r)|.  (11) 

V2  KRm  o<r<t 

Moreover,  if  the  target  normal  acceleration  vanishes,  the  origin  of  the 
closed-loop  system  will  be  exponentially  stable. 

Proof.  The  derivative  of  V  =  \ <o 2  along  the  trajectories  of  system 
(3)  is  given  by 

V  =  coq^-2jCoq+^  -cos(q-<oM)yj  •  (12) 


Theorem  2  shows  that,  with  the  guidance  law  (10),  the  LOS 
rate  coq  converges  to  a  small  neighborhood  of  zero  by  adjusting 
coefficients  f  and  K,  even  though  the  unknown  target  normal 
acceleration  exists.  As  R  becomes  small,  some  smaller  |  or  larger 
K  is  required  to  keep  \coq\  in  a  desired  domain  as  implied  by 
Eq.  (15),  i.e.,  a  large  control  energy  is  needed  to  suppress  the 
target  maneuvers,  since  R  =  0  is  a  singular  point  of  system  (3). 
However,  due  to  the  finite  size  of  missiles  and  targets,  a  successful 
interception  can  be  achieved  as  long  as  the  missile  can  supply 
enough  control  energy  until  R  decreases  to  a  particular  intercept 
value  in  the  whole  process  of  homing  guidance  (Shtessel  et  al„ 
2009;  Shtessel  &  Tournes,  2009). 

3.3.  ICC  law  design 

Since  the  angle  of  attack  must  be  kept  in  a  reasonable  domain 
during  the  end  game,  the  command  profile  a*  can  be  chosen  as  a 
predetermined  bounded  smooth  function.  Taking  the  derivative  of 
ah  (see  Appendix)  yields 

■*  2aTRcoq  raTq 

a*M  =  c - — - 1 - - -  (18) 

cos(q  -  cpM)  R  cos (q  -  tpM) 

where  r  =  R(K  +  —  2VR,  aTR  is  the  projection  of  target  accel¬ 

eration  along  the  LOS,  and  q  is  given  by 
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aM  \ 


vM 


'  VMrc 


'  VMTt 


'  vM’J 


iV*., 


COS (q  -  I Pm)Vm 


aTq 


—  2 

R  Wq 

R 

c0x  + 

fO 

I  1 

0  O' 

■  0 

0 

Aa  0 

Aa  0 

Agc 

Ag, 

0 

i  0 

0  1 

x  + 

0 

0 

_Bl 

*  0 

B,„  0_ 

JSc 

Bgt 

"  0 

0  " 

"  0  " 

— A$c 

— Asc 

Aa 

V 

0 

0 

+ 

0 

_  % 

H  _ 

_ab_ 

where  Co  =  [0  —  1  1  0],  By  choosing 

GMHil 


-AgcBsc  +AscBsc  AstBsc  +AscBst  - 


~k&B& 


x  =  Ax  +  Bu  +  FA 

where  the  matrices  A,  B,  F,  and  A  are  given  by 


'0 

1  0 

O' 

-  0 

0  ' 

A  = 

Aa 

0 

0  Aa 

0  0 

0 

,  B  = 

-A* 

0 

—Age 

0 

_Ba 

0  Bm 

0_ 

_  Bgc 

Bgc  _ 

B  I 

'o 

1  0  ol 

f 

.  [Aa] 

f  =  l 

0 

0  0  lj 

1  ' 

H4 

r,  the  change  of  variables 


T)=X~X 

brings  Eqs.  (24)  into  the  form 


(19) 


H - - - -  I  — —o)q  -  cos(q  -  <pM)  — 

cos (q  -cpM)  \  R  R  / 

where  all  the  state  variables  including  missile  acceleration  are 
assumed  to  be  measurable.  Thus,  according  to  Eqs.  (2)  and  (4d),  we 
define  pitch  rate  command  co*  and  its  derivative  command  co*  as 

=  **  =  i*  +  i-  (20) 

It  can  be  seen  that  if  command  profiles  a*  and  of  are  enforced  by 
control  surfaces,  the  LOS  rate  will  be  1SS  against  the  target  maneu¬ 
vers  according  to  Theorem  2.  Therefore,  we  will  design  a  control 
law  to  drive  a  and  co  to  track  their  commands. 

From  Eqs.  (4),  we  have 

d  =  Aaa  +  Aaco  I  Ajc<5c  1  —  —  Age  8)  -\-  A  a  (21a) 

®  =  Baa  +  Bmco  +  BScSc  +  BStSt  +  Bsc8cc  +  Bsc8ct  +  AB  (21b) 

where Aa  =  ma  +  ^  Am  =  ma  -  ^ ,ASc  =  mSc  +  ^  + 

lSc  _  jsc_  A  r  _  kc  A  _  m  I  ,a^t  I  _  _W_  A  r  _ 

V.,r,  VM’A&c  —  Vurr’Ast  —  +  v2  +  Vu’Agt  ~ 


cos (q  -  (Pm)Vm 


C0B+  -z-B 


(26a) 


»)  =  Ar)  +  Bu  +  Ax*  -  x*  +  FA.  (26b) 

It  is  not  difficult  to  see  that  the  system  can  be  represented  in  the 
form  of  (6)-(7)  with  Hi  defined  by 

cos (q-(pM)a*M  .  v-i 


.  cos (q-(pM)VM 


(27) 


v^,Ba  =  ma+mwma- ma^,Bm  =  mv+ma+ml ,  BSc  =  mSc  + 
^Ico^Sc  Hla  Xc  >  Bgc  —  Tc  )  Bgt  —  lflgt  A  Cl.,, nig,  TTla 

~~~ ,  Ugc  =  •  ,  Aa  =  Aw —  Aa-\~ Aa,  AB  =  TYla  Aa-\-ifiw AM-\- AM. 

Denotingx  =  [a,  a,  co,  cw]T,  on  account  of  Eqs.  (2),  (3)  and  (21),  we 
obtain 


Eq.  (22)  can  be  trans- 


-yz 

\Aaoc*  +Aaoo*  -  a*', 

y^lB^+B^-cb*  J- 

H2  defined  by 

|  i)  =  At]  +  Bu  +  Fv2  +  Fyi  f2g, 

\yz  =  ~c0ri  (  1 

and 

Vj  =  arq,  v2  =  A. 

As  a  result  of  Theorem  2,  system  H  i  is  ISS  against  V\  andy2,  and  the 
following  inequality  holds, 

H(t)l  <  ^ 1  ~  e~2Kt psuP  l^iWI 

+  V 1  -  e~2Kt  sup  |y2(t)|  +  e  Kt|a)?(0)| 

V2 KRm  o  <r<t 

A  a0  ^sup  |ui(r)|j  +  ao  ^sup  |y2(r)|j 

+  A)(K(0)|,t).  (29) 

Hence,  according  to  Proposition  3.1  of  Jiang  et  al.  (1994),  for  the 
output  function  y],  the  inequality 

llyill  <  Y\  (sup  |«i(t)M 

+  ft(l«,(0)|,t)  +  dj  (30) 

holds  for  a  pair  of  class  X  functions  (y",  yf ),  a  class  XX.  function 
Pi,  and  a  nonnegative  constant  dj.  The  main  results  are  stated  in 
the  following  theorem. 

Theorem  3.  Consider  the  interconnected  system  (27)  and  (28).  Take 
control  law  (31)  as  shown  in  Box  I  for  a  >  0,  b  >  0,  and  e  is 
a  positive  constant  to  be  specified.  If  ||A||  is  bounded,  there  exists 
e*  >  0  such  that  system  (27)  is  ISpS  with  respect  to  inputs  Vi  and 
v2for  0  <  e  <  e*. 

Proof.  Substituting  control  law  (31)  and  the  new  variables 


Yi  ^sup  |y2(r)|  j 


(24a) 

(24b) 


»?i 

m  =  — ,  m  =  m,  m  =  - 

into  system  (28)  yields 


h  4  =  n  4 


(32) 


r j  =  -G0rj  +  F  v2  +Fyi 

e 

y2  =  -c0D(e)r) 

where  rj  =  \rji,rj2,rj3,rj4]T,D(e)  =  diag[e,  1,  < 


Go  — 


-a  -b  0  0 

0  0  0  1 
0  0  -a  -b 
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-  Age  Ba  + 

(a*  +  ?)  Bst 

!% 

AmBsc  +  (b^  +  A5c 

lAsc 

u  = 

AgcBg 

AgcBa  + 

c- AgcBgc 
(A*  +  yi)  Bsc 

AgcBgc  —  AgcBgc 

|  % 

AgcBgc  —AscBsc 

Aa>Bgc  +  ^Bm  +  Age 

AgcBgc  —  AgcBgc 

lA% 

r)  —  Pri- 

(31) 

AgcBg 

c  -AgcBgc 

AgcBgc  -AgcBgc 

AgfBgc  —  AgcBgc 

AgcBgc  —  AgcBgc  _ 

Box  I. 


is  a  Hurwitz  matrix.  The  derivative  of  Lyapunov  function  V(fj)  = 
r^Porj,  where  P0  is  the  positive  definite  solution  of  the  Lyapunov 
equation  P0G0  +  GjP0  =  — L  along  the  trajectories  of  system  (33) 
is  given  by 


V  =  --||>7||2  +  2 rjTP0Fv2  +  2 rjTP0Fy,. 

£ 

Substituting  the  inequalities 
>fPoFu2<i||l7||2+^||F||2||Po|rt|%||2 
J7TP0Fyi  <  |ll^7ll2+  |llF||2||P0||2||yi||2 
into  Eq.  (34).  we  obtain 

V  <  ll'7ll2  +  l|P|l2||Pol|2||o2||2 

+  l|p||2IIPoll2llyill2. 

We  will  always  consider  e  <  |.  Since 

^min(Po)ll»?(t)||2  <  V(rj(t))  <  ^max(Po)ll>7(t)||2 
we  have 

ll>7(OII2  <  A3e-AP||7j(0)||2  +  ^(1  -  e-2'() 

x  ^sup  ||y! (r)||2  +^sup  ||u2(r)||2 

and 


(34) 


(35) 


(36) 

(37) 


(38) 


1*2(01  <  II c0 1|  ||D(e)ij(0ll _ 

<  II c0 1|  ||D(e)||J^(t-e-*J*)  sup  ||u2(r) 


f  llcoll  ||D(e)||^(l-e-2 

P||co||v/^=^7lh(0)|| 

Vi  (osup<tl|o2(r)||)  +  yI  ^ 


Due  to  the  form  of  y2 ,  there  exists  a  constant  £q  >  0  such  that 
Eq.  (9)  holds  for  S/  =  0  if  e  <  e*  =  minji,  ej}.  In  this  case,  due 
to  Theorem  1,  we  have 

11*2(011 


<  A(||z(0)||,t)  +  yvi  ^supJutCr)^ 
+  Yv2  ^supj|u2(r)||^  +  d2 


where  z  =  [coq,  rj]T,  yvl  and  yv2  are  class  X  functions,  pz  is  a  class 
XX  function,  and  d2  is  a  nonnegative  constant.  It  can  be  shown 
from  Eqs.  (29)  and  (40)  that2 

KCOI  <  Y'v i  ^sup  |ui(r)M  +  y'v2  ^sup  ||u2(r)||j 

+  ^z(l|z(0)||,  0  +  dm  (41) 

holds  for  a  pair  of  X  functions  ()/-, ,  y’v2),  a  class  XX  function  f)'z, 
and  a  nonnegative  constant  dm.  Therefore,  by  taking  control  law 
(31)  with  e  <  £*,  system  (27)  is  ISpS  with  respect  to  v\  and  v2.  □ 


Substituting  Eq.  (30)  into  Eq.  (39)  yields 

1*2(01  <  y2  ^sup^  ||u2(r)||  j  +  yy2  (yf  ^sup^  |U](T)|  j  +  y{ 

x  ^sup  |*2(r)|^  +  ^i(|cu9(0)|,  0  +  dj] 

+  >02(11^0)11,  0 

<  yl  o  (Id  +  p,)  o  y{  ^sup  |*2(r)|j  +  y*  o  (Id  +  pj"1) 
x  (y”  ^supf  K(r)|)  +  px(\a)q(0)\,  t )  +  df) 

+Y2  ^supj|u2(r)||^  +i62(||J7(0)||,0-  (42) 

A  fact  to  be  noticed  is  that  Eq.  (9)  implies  the  existence  of  a 
nonnegative  real  number  d0  such  that 

{  yl  o  (Id  +  pi)  o  Yj  (s)  <  (Id  +  p2yx  (s)  +  do  y 
l y{  o  (Id  +  p2)  o  y2y(s)  <  (Id  +  px)-' (s)  +  do  ’ 

with  d0  =  0  if  S|  =  0.  Thus,  Eq.  (43)  holds 
1*2(01  <  (Id  +  ft)-1  ^sup  |*2(r)|J  +  d0  +  y2  o  (Id  +  pj-1) 
x  (y"  ^sup  | j  +  y8i(|wg(0)|,  t)  +  d^j 
+  Y2v  ^sup  ||i>2(r)||^  +  ft(||f/(0)||,  0 

<  (Id  +  p^1)  (ft  (||  >7(0)  II,  0  +  d0  +  yl  o  (Id  +  Pxx) 
x  (V;  ^sup^  |ut(r)|^  +  Px(\(Dq(0)\,  0  +  dfj 

+  y2v  ^sup  ||u2(r)||^.  (43) 


2  For  any  class  X  function  y,  any  class  Xoa  function  p  such  that  p  —  Id  is  of  class 
Xoo,  and  any  nonnegative  real  numbers  a  and  b  we  have 

Y(a  +  b )  <  y(p(a))  +  y(p  °  (P  -  *0_10>)). 
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(a)  Relative  range. 


(b)  LOS  rate. 


(c)  Angle  of  attack. 


(d)  Pitch  angle. 


Fig.  2.  The  performance  of  the  proposed  ICC  law  in  case  1. 


Substituting  Eq.  (43)  into  Eq.  (29)  yields 

H(t)|  <  ao^sup  |t>i(r)|^  +  a0  ^(/d  +  p2-1) 

x  ^(1117(0)11,0+ ffg +  }<*  ^supt||t;2(T)||j 

+  y2  o  (Id  +  pj-1)  | '  y tv  ^sup^  |ui(t)|^ 

+  fh(lcog(0)l,  t)  +  +  A,(|<a,(0)|,  t).  (44) 

Since  d0  =  0  holds  for  0  <  e  <  e*,  and  y2”,  y2  -#>  0  as  s  — >  0,  it 
can  be  seen  that  the  right-hand  side  of  (44)  approaches 

a0  (osupf  +  «0  °  Qd + 

+A)(H(0)l,0- 

as  e  -f  0  for  bounded  V\,  v2  and  d1,  which  shows  that  for 
sufficiently  small  e  the  influence  of  u2  and  dj  on  coq  will  be  close 
to  zero.  Besides  that,  since  a0(x )  =  Vf  —  e~2Ktx,  iq  can  be 
suppressed  by  adjusting  f  and  K. 

According  to  Theorem  3  and  Eq.  (23),  the  ICC  law  can  be  written 
as 


Remark  2.  The  form  of  term  Prj  is  similar  to  a  PD  controller,  which 
aims  at  providing  a  sufficiently  fast  tracking  of  the  command 
profiles  to  guarantee  stability  of  the  integrated  system. 

4.  Simulation  results 

This  section  presents  simulation  results  for  the  closed-loop 
system  with  the  proposed  IGC  law  (45).  In  this  simulation  study 
the  constant  missile  speed  is  assumed  to  be  VM  =  380  m/s2.  The 
missile  model  parameters  are  la  =  800  m/s2,  lgc  =  h,  =  60  m/s2, 
ma  =  -100  s  ~2,mm  =  -5s_1,mac  =  50  s  ~2,mSt  =  -50  s-2, 
and  rc  =  rt  =  0.03  s-1.  The  initial  missile  attitude  and  control 


fins  are  a(0)  =  0.05  rad,  if  (0)  =  0.8  rad,  eo(0)  =  0,  and  <$c(0)  = 
<$t(0)  =  0.  The  actuator  commands  are  assumed  to  be  limited 
by  |5C|  <  |,  | (5, |  <  |,  and  the  maximum  missile  acceleration 
is  =  100  m/s2.  The  initial  variables  of  relative  motion  are 
R( 0)  =  3000  m  and  q( 0)  =  The  initial  target  speed  and  fight 
path  angle  are  VT  =  200  m/s  and  <pT(0)  =  1.32  rad.  It  is  easy  to 
calculate  the  initial  LOS  rate,  ojq(0)  =  0.03  rad/s.  Here,  we  define 
the  command  profile  a*  to  be  a  hyperbolic  tangent  function 
a*  =  0.04tanh(a/0.04)  =  0.04^^=^^,  which  is  a  smooth 
function  bounded  by  0.04  rad.  The  LOS  rate  is  taken  as  a  first- 
order  lag  system  with  a  time  constant  0.005  s  with  the  additive 
Gaussian  noise  with  mean  zero  and  standard  deviation  0.1  mrad, 
and  other  feedback  states  are  taken  with  the  same  noise  as  well. 
Uncertainties  Aa  and  A,„  are  set  as  Aa  =  0.03sin(t)  rad/s  + 
Ara,  Am  =  0.3  cos(t)  rad/s2  +  A'm,  where  Ara  and  Arm  are  random 
uncertainties  of  Gaussian  noises  with  deviations  of  0.01  rad/s 
and  0.2  rad/s2,  respectively.  The  contribution  of  gravity  to  lift  of 
gcos(i?-a)  js  a|SQ  ta|<en  jnt0  consideration  in  the  simulation.  Because 

vM 

the  angles  using  ‘radian’  may  seem  to  be  small  but  the  actual  angle 
may  be  large,  the  unit  of ‘degree’  is  used  in  all  the  simulation  figures 
related  to  angles. 

Case  1:  Suppose  that  the  target  does  not  maneuver.  The 
guidance  and  control  coefficients  are  chosen  as  K  =  5,f  =  l,e  = 
0.01,  and  a  =  b  =  1.  The  simulation  results  are  shown  in  Fig.  2. 
Fig.  2(a)  and  (b)  give  curves  of  the  target-to-missile  range  and  LOS 
rate.  When  homing  guidance  begins,  coQ  rapidly  converges  to  a 
small  neighborhood  of  zero.  The  miss  distance  is  only  about  R  = 
0.05  m,  which  means  a  successful  interception.  The  response  of 
missile  dynamics  during  the  process  of  homing  guidance  is  shown 
in  Fig.  2(c)-(f).  At  the  beginning  of  the  homing  guidance,  the  missile 
acceleration  is  relative  large  so  as  to  adjust  flight  direction  to  nullify 
LOS  rate,  and  when  LOS  rate  converges  to  a  very  small  value,  missile 
acceleration  together  with  angle  of  attack  and  pitch  rate  is  kept 
around  zero.  Fig.  2(g)  and  (h)  illustrate  canard  and  tail  deflections. 

Case  2:  Assume  that  the  target  escapes  with  acceleration  of 
aTq  =  60cos(|t)  m/s2.  The  guidance  and  control  coefficients  are 
chosen  as  K  =  f,  f  =  0.1,  e  =  0.01,  and  a  =  b  =  1.  The  sim¬ 
ulation  results  are  shown  in  Fig.  3.  From  Fig.  3(a)  and  (b)  we  can 
see  that,  when  homing  guidance  begins,  coq  rapidly  converges  to 
a  small  neighborhood  of  zero,  and  can  be  kept  near  zero  although 
the  target  normal  acceleration  always  exists,  until  the  range  R  de¬ 
creases  to  an  extremely  small  value  (R  ~  0.05  m).  Compared 
with  Fig.  2(b),  the  deviation  of  LOS  rate  from  zero  is  a  little  larger 
due  to  target  maneuvers,  but  a  successful  interception  can  also  be 


